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SUMMARY 


This report describes an experiment designed to determine the dynamic behavior of 
a drag -force anemometer in high-frequency, unsteady flow. In steady flow the output of 
the anemometer is proportional to stream velocity head and flow angle. Fluid mechan- 
ics suggests that, in unsteady flow, the output would also be proportional to the rate of 
change of fluid velocity. It was determined from the experiment described herein that 
effects due to the rate of change of fluid velocity are negligible for the probe geometry 
and frequencies involved. 


INTRODUCTION 

The drag-force anemometer considered in this report is a cantilevered beam with 
attached strain gages . It is capable of measuring both static and dynamic velocity head 
and flow angle. 

In contrast to the hot-wire anemometer, which is an energy -balance device, the 
drag-force anemometer has an output that is proportional to the force exerted on the 
beam by the moving fluid. For steady-state conditions the force is just that due to the 
velocity head of the fluid. However, fluid mechanics suggests that, if the velocity is 
not constant, an additional term is present that is proportional to the rate of change of 
fluid velocity. 

If this term were not negligible, there would be a spurious response at high fre- 
quencies since the velocity derivative would be higher there. It would also be more 
difficult to extend the frequency response of the beam electronically (ref. 1). 

This report describes the effect of a velocity derivative term and the experiments 
performed to detect that effect. 
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frequency, Hz 

natural frequency of drag -force anemometer, Hz 

I0 (h/10) 

spacing by which equal -error curve misses specified data points 
offset parameter in curve -fitting scheme 
output amplitude of drag-force anemometer 
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value of 1 data point 
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velocity of fluid, m/sec 
steady part of flow, m/sec 

electrical output of drag-force anemometer, mV 
dc part of signal, mV 
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curve-fitting parameter, 1 - 2£ 
velocity derivative parameter, (C/2 k) 
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damping coefficient of drag-force anemometer 

o 

velocity head coefficient, mV/(m/sec) 
turbulence amplitude of moving fluid 
density, kg/m 

phase shift in drag-force anemometer output 
27rf, sec ^ 

27rf n , sec * 


THEORY 

The drag-force anemometer considered in this report is a cantilevered beam with 
strain gages attached near the fixed end (fig. 1) . It is a sdfcond-order system (ref. 2) 
and is described by the second -order differential equation 
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The right member of equation (1), Dj-j(t), denotes the drag on the beam caused by the 
moving fluid. The natural frequency f = co n /27r is a function of the length, thickness, 
and composition of the beam; the damping coefficient £ is a function of the composition 
only; x(t) is the output signal. 

For steady flow the drag is proportional to the exposed area of the beam and 
to the velocity head l/2pU 2 . That is, for a steady flow of velocity 

D D ( t )=C D A D (l/2puj) (2) 

The constant of proportionality is called the drag coefficient. 

The theory of fluid mechanics suggests that for unsteady flow, where U is not 
constant, another term is present that is proportional to the rate of change of the fluid 
velocity (ref. 3). For this case equation (1) is rewritten as 
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Written in this way, both k and C have the same units, mV/(m/see) . The case 
where C = 0 corresponds to the absence of the velocity derivative term. 

The theoretical response of the probe can be examined by assuming that the flow 
velocity U(t) is given by 


U(t) = U Q (1 + X sin cot) for X « 1 (4) 

which is a steady flow of velocity plus a small oscillatory turbulence of magnitude 
\U Q and frequency co/27r. 

With U(t) given by equation (4) 

U 2 (t) = U 2 (l + 2X sin cot + X 2 sin 2 cot) 

=* U 2 (l + 2X sin cot) 
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to first order in X, and dU/dt = AU Q co cos cot. The expression for the beam output 
(eq. (3)) becomes 


-iL + — x + x = /cuf + 2XkU 2 sin cot + \CU^— cos cot 
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The solution of equation (5) can be written in the form 
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is the amplitude of the output as a function of frequency and 
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( 8 ) 


is the phase shift. The dimensionless parameter /3, which appears in equations (I) 
and (8), is defined as 


1 


(9) 


It indicates the relative strength of the velocity derivative effect. 
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As can be seen from equation (7), the output of the beam depends on the frequency as 
well as on the parameters /3 and £. Response curves for various values of /3 are shown 
in figure 2 for the highly underdamped case where £ = 0.01. The quantity plotted is 


Relative output (dB) = 20 log 

x(0) 


= 10 log 



( 10 ) 


EXPERIMENTAL SETUP 

One of the drag-force anemometers that was tested is shown in figure 1. The 
probe consists of a commercially available silicon beam 0.25 millimeter thick with an 
unsupported length of 2. 5 millimeters, the outer 1.5 millimeters of which is exposed 
to the flow. A four-arm, diffusion-bonded strain-gage bridge is on one side of the 
beam at the base. Because silicon is brittle like glass, the anemometer is fragile and 
care is required in handling it. The natural frequency of the beam in the first bending 
mode is about 40 kilohertz. Details of the anemometer are given in reference 2. 

Flow excitation to the anemometer was supplied by a 9.5-millimeter-radius free 
jet that was operated at a Mach number of 0. 3. The jet turbulence intensity was about 
2 percent. The silicon beam was mounted 5 millimeters from the end of the free-jet 
nozzle. To determine the true frequency spectrum of the free jet, a hot-wire anemom- 
eter with a frequency response of at least 80 kilohertz was used . The hot wire was 
placed at the same location in the free jet as the beam, with the wire parallel to the 
plane of the jet exit. 

The output of the anemometer was fed into a spectrum analyzer that had a band- 
width of 20 hertz to 50 kilohertz. The analog output of the analyzer was connected to 
the plotter to produce hard copies of anemometer output as a function of frequency. 


RESULTS 

Figure 3 shows the output signal spectrum from the drag-force anemometer, and 
figure 4 shows the spectrum of the open jet obtained with the hot-wire anemometer. 
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Since the hot-wire anemometer has a frequency response of at least 80 kilohertz, it is 
assumed that figure 4 is also the true spectrum of the jet. Since the true spectrum of 
the jet is not flat with frequency, an appropriate correction must be made to the output 
signal spectrum of figure 3 in order to obtain the desired frequency response of the 
drag -force anemometer. The amount of the correction at each frequency is equal to the 
number of decibels that must be added to the jet spectrum (fig. 4) in order to bring the 
curve up to zero. Figure 5 shows the corrected spectrum for the drag -force anemom- 
eter, which is identical to its frequency response curve. 

The theoretical response equation (eq, (10)) was then fitted to the curve of figure 5 
by using the "min-max" method in order to determine the parameters /3 and £. The 
details of this procedure are given in the appendix. In using this procedure the natural 
frequency f is considered to be known since for a highly underdamped system the 
peak occurs at nearly f (ref. 2) . 

The experimental and fitted curves are shown in figure 6. The maximum deviation 
h between the experimental and analytical curves was 0.9 decibel. The parameter h 
is the figure of merit for the min-max method. 

The fitted curve corresponded to £ and /3 values of 0 . 03 and 0 . 2, respectively. 

A second drag -force anemometer was tested and gave values for £ and (3 of 0.02 and 
-0.3, respectively. The figure of merit h for this second curve fit was 1.9 decibels. 
Of course, /3 cannot really be negative since it is a squared quantity, but the indication 
is that /3 « 1 for both anemometers tested. Therefore the velocity derivative effect 
can be ignored, at least up to f/f = 1. 25, or to about 50 kilohertz for the high- 
frequency silicon beams considered herein (f =* 42 kHz) . 


CONCLUDING REMARKS 

The results of theory and experiment on miniature, beam drag-force anemometers 
show that the velocity derivative term in the drag-force function can be neglected for a 
wide range of frequencies. For the case shown, where silicon beams with a natural 
frequency f of 42 kilohertz were used, this velocity derivative term could be neglec- 
ted to at least 50 kilohertz. This range should certainly be adequate since in most 
cases the user is interested only in frequencies less than about 0.25 f 

Electronically extending the frequency range of the anemometer should present no 
problem . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 20, 1980, 

505-04. 
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APPENDIX - CURVE FITTING 


The curve to be fitted has the form 


1 + p(t/f 

R = 10 log V_JV + k (Al) 

[( f/f n) 2 - if + « 2 ( f / f n) 2 

The argument of the logarithm is from equation (7); the logarithm occurs because the 
spectrum analyzer output is in decibels. The extra term K is added because the zero 
point is arbitrary. 

The usual method of curve fitting is least squares, where the fitted curve satisfies the 
criterion that the sum of the squares of the deviations between the analytical and experi- 
mental curve is a minimum. Because of the nonlinearity of equation (Al), it was neces- 
sary to iterate to find (3 and £ by using the least-squares method . However, the itera- 
tion diverged. Furthermore it tended to produce a good fit only near the peak at f/f =1. 

For these reasons the least-squares approach was discarded. Instead a curve- 
fitting technique called the min-max method (ref. 4) was used. In this method the curve 
fitted to the data satisfies the criterion that the maximum error between the two curves 
be a minimum. The existence and uniqueness of the min-max curve are discussed in 
the reference, as are the equal-error property and the exchange algorithm. 

The equal-error property is that property of a min-max curve by which it deviates 
from the data with alternating sign by an amount h, the maximum error, exactly n + 1 
times, where n is the number of parameters that describe the curve. The exchange 
algorithm makes use of this property in actually fitting the curve to the data. The pro- 
cedure starts with the fitting of an equal-error curve to n + 1 arbitrarily selected data 
points. If no other data point has an error larger than these n + 1 points, the desired 
curve has been found. If one or more other points do have larger errors, the point hav- 
ing the largest error is used to replace one of the original points in such a way that the 
errors are still of opposite sign. A new equal-error curve is then fitted to this new set 
of n + 1 points, and the process continues until no data point has an error greater than 
the n + 1 points used in constructing the curve . 

In this case, the function to be used to fit the curve of figure 5 is given by equa- 
tion (Al), so that there are three parameters: £, /?, and K. At this point it is con- 
venient to make the change of variables 


7 


L 





z = 


n; 


(A2) 


OL - 1 - 2q J 


In terms of these new variables, equation (Al) becomes 


R(z) = 10 log( . 1 + ) + K 


(A3) 


yz - 2az + 1, 


Now let z. be the locations of four arbitrary data points, and r^ their magni- 
tudes, for i = 1, 2, 3, 4. The equal-error curve is to miss each of these points by an 
amount h, with alternating sign. That is, R(z i ) - = ±h, so that there are four equa- 

tions in the four unknowns a, p, K, and h. These equations are 
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These equations can be simplified by dividing through by 10 and then taking the antilog 
of each equation. For instance, the first equation becomes 
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in 1 * iptK/lO) = 10 (r l / 10 ) 10 (h/10) 

1 - 20^ + Z 2 


From the definitions 


= io (k/io> 


s. = 10 


(r./10) 


H = io (b/ 10 ) 


this equation becomes 


y + /3' z i 


1 - 20^ + z^ 


where /?’ = y|3. Likewise equations (A4) simplify to 


y + z + 2Hs^z^c>' - SjH^l + z^ 
Hy + Hz^’ + 2s 2 z 2a = s 2 (l + z|) 
y + Zg/3 1 + 2Hs 3 z 3 a = SgH^l + z^ 
Hy + Hz 4 /3» + 2s 3 z 4 q; = s 4 jl + z 


If we call D the determinant of the coefficients, equations (A 6) have the solution 
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and D is 


D = H 2 [s iZi ( Z3 - z 2 ) + s 3 z 3 (z 2 - zj] - s 2 z 2 (z 3 - Zl ) 

These quantities can be substituted back into equations ( A6 ) to find H. The result is 

2 

the quadratic equation for H 
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+ Sl s 2 (l - Zl z 2 )] + (z 4 - z t ) (z 3 - z 2 ) [s 2 s 3 (1 - Z2 Z 3 ) + V4 (1 - z jZ4 )] | 




2 

Of the two values of H that result, the positive one is chosen. If they are both posi- 
tive, the one that yields the smaller value of h is used. Once H has been found, 
equations (Al) give the parameters a, (3, and y for the four points chosen. The 
process continues as explained until all the data points are encompassed within the 
equal -error curve . Then all the parameters have been found for the curve that fits the 
data with the smallest maximum error, and the value of h found in the process gives 
a measure of how well the data have been fitted . 
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Frequency, kHz 

Figure 3. - Output signal spectrum of typical drag-force 
anemometer in open jet. 



Frequency, kHz 

Figure 4. - Output signal spectrum of hot-wire anemom- 
eter in open jet. 




Frequency, kHz 

Figure 6. - Experimental and fitted frequency response 
curves. 
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